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Abstract 

We propose an energetic interpretation of stochastic processes described by Langevin 
equations with non-uniform temperature. In order to avoid Ito-Stratonovich dilemma, 
we start with a Kramers equation, and derive a Fokker-Plank equation by the renor- 
malization group method. We give a proper definition of heat for the system. Based 
on our formulations, we analyze two examples, the Thomson effect and a Brownian 
motor which realizes the Carnot efficiency. 
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1 Introduction 



Recently, Brownian motors have attracted considerable attention stimulated 
by research on molecular motors. A Brownian motor, which appeared in Feyn- 
man's famous textbook for the first time as a thermal ratchet [1], is the ma- 
chine which can rectify thermal fluctuations and product a directed current. 
Until now, several models of a Brownian motor have been proposed. The 
most received and extensively researched system consists of Brownian parti- 
cles subjected to a Gaussian white noise whose intensity is characterized by a 
temperature for the environment [2] . There are other types of Brownian motor 
models. In particular, it is known that a Brownian particle under non-uniform 
temperature profile moves directedly on the average when the potential profile 
synchronizes the temperature profile but they are not same shapes [5]. Here, 
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the motion of the particle was assumed to be described by Langevin dynamics 
supplemented with a state dependent noise. 

Until now, the energetic consideration for a Langevin dynamics with uniform 
temperature has been developed [3,4]. In fact, the efficiency for the Feynman's 
ratchet was discussed based on the notions of heat and work which are defined 
properly in a consistent way with thermodynamic laws. Then, the question 
we address here is how to define heat and work in non-uniform temperature 
systems. Thus, our aim of this paper is to extend the energetic consideration so 
as to apply to such systems. It will turn out that this question is not trivially 
solved owing to the nature of the state dependent noise. 

In this paper, so as to avoid the complications for the state dependent noise, 
we first derive a Fokker-Plank equation from a uniquely defined Kramers equa- 
tion by employing the perturbative expansion. There are many perturbation 
methods to derive a Fokker-Plank equation [6,7]. Since we treat an delicate 
problem, the most formal theory without heuristicity might be best. Thus, we 
adopt renormalization group (RG) method, whose applicability to differential 
equations has been established recently [8]. Then, following the Sekimoto's 
formulation [3], we make an energetic interpretation of the system. After that, 
as applications, two kinds of non-uniform temperature systems are analyzed. 
As the first example, we study the Thomson effect for non-interacting parti- 
cles and show that the Thomson coefficient is 1/2. As the second, we calculate 
the work efficiency for the Biittiker's model and show that the Carnot limit is 
realized under the quasistatic condition. 

This article is organized as follows. In section2, we review a sort of dilemma 
caused by the nature of a state dependent noise. In section3, the Fokker-Plank 
equation is derived from the Kramers equation by employing the renormaliza- 
tion group method. In section4, we define heat for the system described by the 
Fokker-Plank equation. In section5, we analyze above two examples. Section6 
is devoted to a summary and conclusions of this work. 



2 Ito-Stratonovich dilemma 

In this section, we review a dilemma caused by ill-defined nature of a Langevin 
equation with a position dependent noise intensity [9] 




(1) 



mm) = - n , 
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where 7 is a friction constant of a Brownian particle, / is an external force, 
T(x) is a temperature profile and we have assumed that a Boltzmann con- 
stant k is equal to the unity. Owing to space dependence of a noise intensity, 
the following difficulty arises. Because the noise has no correlation time, it is 
not yet clear which x value one has to use in the function T(x) 1/2 in Eq.(l). 
If, for instance, £ is considered as a sum of peaked functions with no width, 
the stochastic variable x will jump at every time when such a peaked func- 
tion occurs. The question then arises: which x value must use in T(x) 1 / 2 ? 
One may use the value x just before the jump, after the jump, or some value 
between these two values. From a purely mathematical point one cannot an- 
swer this question, but one has to use a prescription as additional specifica- 
tion. Among some prescriptions, there are two representatives called Ito and 
Stratonovich. Under the Ito prescription, the value of x is inserted before the 
jump, while the mean value between before and after the jump is opted for, 
under the Stratonovich prescription. The Fokker-Plank formalism explicitly 
shows this prescription dependence as follows: A Fokker-Plank equation with 
the Stratonovich prescription is written as 

9(j) ld (T 9 f^ l9T U ™ 



while a Fokker-Plank equation with the Ito prescription is written as 

^ = IA( T A_ / + ^U (4) 



Such a dependence of an apparently unique equation on prescriptions is called 
"Ito-Stratonovich dilemma", and their difference 

itekte)* (5) 

is called a spurious flow. Owing to the existence of this dilemma, we cannot 
proceed to discuss energetics, because descriptions of energetics are subject to 
prescriptions. 

However, it is possible to avoid this dilemma by considering the effect of 
inertia of a Brownian particle. An underdamped type of a Langevin equation 
is written as 



x = v , 

(6) 
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where we have assumed the mass of a Brownian particle m is equal to the 
unity. From Eq.(6), we see that each delta peak in £(t) causes the jump in 
v(t), but no jump in x(t) and, therefore, T(rr) 1 / 2 . Since does not cause a 
jump in the noise intensity, the solution is uniquely determined for the given 
£(£). As a result, the evolution of a distribution function P is uniquely denoted 
by 

which is called a Kramers equation. Since there is no complication for the 
prescription dependence, we discuss Eqs.(6) and (8). 



3 Perturbation and Renormalization 

In this section, employing a renormalization group method, we derive a Fokker- 
Plank equation from a Kramers equation with non-uniform temperature. Now 
we restrict ourselves to the one dimensional case so as to simplify the analysis. 
First, we rewrite Eq.(8) as 




We consider a case that 7 is sufficiently large, that is, a constant I/7 is re- 
garded as a small parameter. To emphasize smallness of the constant, we 
rewrite I/7 as e. Also, for the later convenience, we rescale the time variable 
t to et. The resultant equation takes the form 

This rewriting of the time variable is not essential to the problem, but neces- 
sary for us to employ the renormalization group method properly. 

Now we seek solutions under the condition e — > 0. We first expand the distri- 
bution function with e in such a way that 

P(x,v,t) ~ P {0) (x,v,t) + eP {1) (x,v,t) + ■ ■ ■ . (11) 
Identifying terms of the same order, we obtain 
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O(£ « ): | ; (.„ + r(a;) |.)p(»._|p(»» = , (12) 



We may successively solve these equations from lower order. As preparation 
for solving them, we define the Fokker-Plank operator as 

£ " S |(« +T |)- < 15 > 



and introduce the eigenfunctions of the operator, which are written in the 
form 

L " {x ' v) = ik) Hn{ 7m )eM ~^ (16) 



where H n is the n-th Hermite polynomial. The characteristic equation is given 
by 

L>FpL n (x,v) — —nL n (x,v). (17) 
We also expand p( m ) in terms of eigenfunctions: 

oo 

P {m \x,v,t) = Y,a ( r ) (x,t)L n (x,v). (18) 

n=0 

Now, let us solve Eq.(12). Decomposing left-handed-side of Eq.(12) in terms 
of the eigenfunctions, we obtain 

^ag\x,t) = -nag\ (19) 

We notice that a n , (n ^ 0), decays exponentially in time. Thus, discarding 
such terms, we obtain 

P(0, = # eXp (-^' (20) 
where <f>(x) is an arbitrary function in x. 
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The physical picture behind this equation is as follows : In the fast time scale 
t ~ 1, the velocity distribution relaxes to the Maxwellian rapidly and space 
distribution stays just as it is. 

We next proceed to calculate the higher order terms. Substituting the solution 
Eq.(20) into Eq.(13), we obtain 



T'(j) v 3 



dt 



2T 2 y^T 



v 2 (86 
^dx 



f_± 
T 



8T 



T'6L 3 + 
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T'0x v 
2T J 
V 



2T 



T + T 



(21) 
(22) 



where T' denotes dT/dx. Eq.(22) leads to a set of equations 



3a?> - la® 
3 dt 3 



8T 



n {1) 
a i 



d 
d 



2T/90 
2 



/0 r 

T T 



-na^--a^ = (n^l,3). 



(23) 

(24) 
(25) 



Neglecting terms decaying exponentially in time, we obtain 



3 8T ^ 3 2 t T ' 



(26) 



= e ~2T _ (_L - ^21 _|_ 0) e -2T. (27) 



6T 2 v^ V&r T IT^'^jT 

Then we substitute into the 0(e 2 ) equation so that the right-handed- side 
of the equation takes the form 



>■ 8x ' dv> >\d x T IT > yff 



2 

e - 2T 



d (86 f6 T' \ v 2 ^ 

1 — — — -I 6) e~2T 

dx^dx T 2T V ) 

T' (86 f6 T' \ v 2 „2 
H \-zr- — H <M^e~^ (28) 

T'6f v 2 .2 / 50 /0 j" w 2 ^ 2 

— e~zr + — I — — — H <i>l e - ^ 1 

2T 2 ^ ^T^a; T 2T ' ^/T 
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^d_(^±\^_ -|i _ T' ,d<f> _ f$ V s v 4 



T'<j>f v 4 _Z T ,2 v 4 _vL T /2 <p v 6 _vl 

-I P 2T -|~ P 2T — p 2T 

6P 3 v/T 12T 3 vT 12T 4 vT 



In the same way as before, we derive equations for affl. In particular, 
satisfies 

- |4 2) = V>, (29) 



where 



* s /(»s + 'l) pW< "' (30) 



Therefore, P^ is expressed as 

1 t, 2 

p( 2 ) = — e 2 ipt—=e~%r + 0(e) (regular term), (32) 
VP 



and the asymptotic expansion up to the second order is written as 

1 v 2 

P = (<f>(x) - e 2 ipt)—=e~%r + 0(e) (regular term). (33) 
VP 



This shows that P has a secular term which is unbounded in time. Owing to 
the existence of the secular term, the perturbation is valid only in a short time 
interval satisfying e 2 t <^\. Improve the result of naive perturbation, we apply 
the renormalization group method. We introduce a time r when the initial 
condition for is given. That is, we rewrite Eq.(33) as 

1 v 2 

P = (<f)(x, t) — e 2 ip(t — r))^^e~2T + 0(e) (regular term), (34) 

VP 



Since P should not depend on the choice of r, we impose the condition 



which is called a renormalization group equation [8]. This yields 

at 



(36) 



Restoring the time scale to the original one and using Eq.(31), we finally 
obtain 



d(f)(x,t) 1 d 



dt 7 dx v 9a; 



(T(x)-£ - f 4 + ?<!>)= 0. (37) 



4 Energetics 



Equation (37) takes the same form as a Fokker-Plank equation given by Eq.(3). 
Therefore, an overdamped Langevin equation corresponding to Eq.(37) is sim- 
ply given by Eq.(l) when the Ito prescription is assumed. However, we notice 
that this simpleness dose not ensure the uniqueness of prescription. We can 
consider that Eq.(37) corresponds the Langevin equation given by 

T^V^ + Z-^, (38) 



when the Stratonovich prescription is assumed. Although we have two ap- 
parently different forms of Langevin equations, we restrict ourselves to the 
Stratonovich prescription because the Stratonovich prescription is more suit- 
able for description of energetics, as we will see below. 

The revision term which we call the third term of right-handed-side of Eq. (38) 
is relevant when one faces energetics. Following Sekimito [3,10], we define heats 
for the systems described by an overdamped and an underdamped Langevin 
equation with a position dependent noise intensity as the following stochastic 
integrals 

Qovcr = - J{lX- sjT{x)^dx{t), (39) 
Qundcr = ~ J{lV - ^T\x)i}dx{t) , (40) 



where the integral is assumed in the Stratonovich sense. Substituting Eq.(38) 
to Eq.(39), we obtain the decomposition of the heat of the overdamped case 

Qovcr = Qu + Qt, (41) 
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Qu = -J fdx(t), 



2 dx 



dx(t), 



(42) 
(43) 



where Qu is the part which an external force contributes to and Qt is the 
part which position dependent temperature contributes to. On the other hand, 
substituting Eq.(6) to Eq.(40), we get the decomposition of the heat of the 
underdamped case 



Q under — Qll Q Ai 

Q A = J vdx(t). 



(44) 
(45) 



Here, we show that Qt can be derived as the first order contribution to Qa for 
the high friction limit. Taking an ensemble average about both-handed-sides 
of Eq.(45), we get 



(vdx) 



d 



(46) 
(47) 
(48) 



Neglecting anything but the lowest contribution, we obtain 



(Q A )~±Jdt^Jdx J dvv 2 p(°\x,v,t) 



-- I dt I dxT(x)d x J(x,t) 



1 f r dT 1 r r 

= - dt dx—J(x,t) - - dt T(x)J(x,t) 

Zj J J CJX Zi J 

= {Q T )-\jdt[T(x)J{x,t) 



XR 
XL 



XL 



(49) 
(50) 
(51) 
(52) 
(53) 



where the system has been assumed to be an interval [xl,xr\. The surface 
term may be interpreted as the energy flowing from the outside of the system. 
Then, as far as the energy exported form the heat bath is concerned, we have 
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only to pay attention to the bulk term. Thus, we conclude that Qt is the 
lowest order contribution to Qa, that is, 

Qa ~ Qt- (54) 

We shortly account for the meaning of the revision term. We saw in section3 
that the velocity distribution function rapidly relaxes to the position depen- 
dent Maxwellian in the high friction case. The averaged energy of each position 
H is written as 

H(x, t)= J {^v 2 + U(x))P(x, v, t)dv (55) 
~ f(-v 2 + U(x))(f)(x,t)Lo(x,v)dv (7^00) (56) 
= (±T(x) + U(x))<j>(x,t). (57) 

This equation implies that the kinetic energy contributes to the effective po- 
tential via the equipartition law. Thus, the effective potential may lead to the 
revision term. 



5 Examples 

In this section, we investigate two examples as application of the energetic 
interpretation made in section4. 

Ex. 1) Thomson effect 

The Thomson effect, one of the thermoelectric effect, is the phenomenon that 
some heat is generated besides Joule heat when a current flows in a conduc- 
tor at which a temperature gradient is charged. The generated heat is called 
Thomson heat. According to phenomeno logical laws, Thomson heat generated 
per unit time between two points of the conductor, written as Q, takes the 
form 

Q = tJAT, (58) 

where AT is the temperature difference between two points, J is a current 
and r is a material constant called "Thomson coefficient" [11]. 
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Now let us calculate Thomson heat from a kinematic point of view. We assume 
the system in question is an open interval I = (0, 1) on an real axis and 
particle reservoirs are points 0, 1. Next we give a few conditions. We assume 
that the particle number distribution is realized with P , Pi in each reservoirs, 
temperature T(x) is externally given and there is no mechanical potential. 
Using these conditions and Eq.(37), the problem we solve is written as 



P(0) = P P(l) = Pi. (60) 



Now we restrict ourselves to a steady state. From Eq.(59), a steady state 
distribution P s satisfies the relation 



^TP, = -J, (61) 



where J means a particle number flux and is a constant whose value is decided 
by boundary conditions. The solution of Eq.(61) takes the following form 

T(x)P s (x) = -Jx + T(0)P (0<x<l), (62) 
J = T(0)P -T(1)P 1 . (63) 

The heat generated in the system is given by 



Q=\lfjAt). (64) 



We take an ensemble average for both-handed-sides of Eq.(64), and then we 
obtain 



1 r } dT 

(Q) = -Jdtjdx—J, (65) 
o 

= \ AT Jt, (66) 



where AT is defined by 

AT = T(l) — T(0), (67) 
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and we have used the fact that J is a constant. Differentiating both-handed- 
sides of Eq.(66) with t, we finally obtain 



(Q) = \*TJ. 



(68) 



Comparing this relation with the the phenomenological law Eq.(58), we obtain 



Notice that, in our model, interactions between particles are not considered. 
If the effect of the interactions are taken into account, the value of Thomson 
coefficient would be shifted from 1/2. 

Ex. 2) Biittiker's model 

Biittiker derived a rigorous steady solution of Fokker-Plank equations with a 
position dependent temperature [5]. Let us first review his results by employing 
Eq.(37). It is assumed that the system is defined over an interval [0, h] with 
periodic boundary conditions. It is also assumed that the force / is given by 



This model has a steady current solution expressed by 



where Z is a normalization constant [5]. Equation (71) implies that a finite 
current exists when U (x) and T(x) satisfy the condition 



Now, in order to extract the work from this system, a load is superimposed. 
That is, the form of the mechanical potential is assumed as 




(71) 




(72) 



U{x) = U (x) + Lx (L > 0), 



(73) 
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where Uo(x) is a periodic function with a period h, 

U (x + h) = U (x), 



(74) 



and L corresponds to the load per unit length. When the temperaute T and 
a periodic part of the potential U satisfy the inequality 



h 



T(x) dx 
o v ' 



a Brownian particle climbs the slope of potential so that it works against the 
load. For simplicity, we assume the temperature function as 



T(x) 



Ti :xE [0, a] 

[ ' J (76) 

T 2 :xe [a,h], 



where Ti > T 2 . Let us calculate the work efficiency defined by 



where (Qh) is averaged heat extracted from the hotter heat bath whose tem- 
perature is Ti, and (W) is averaged work done by the particle against the 
load. Since the hotter heat bath is located in the interval [0, a], Q H is written 

as 



z(t)€[0,a 



Taking an ensemble average, we get the denominator of Eq.(77) 



' dU 



(Qh}= JdtJ dx—J, (79) 
o 

= {U(a)-U(0)}Jt. (80) 

The work against the load is defined by 

W = Lj dx{t). (81) 
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(For more rigorous definition of work, see [3] and [12].) Then we can calculate 
the numerator Eq.(77) as 



(W)=L(Jdx(t)), (82) 

= LhJt, (83) 
= {U(h) -U(0)}Jt. (84) 

From Eq.(80) and Eq.(84), the efficiency is expressed as 

_ U(h) - U(0) 
1 - U{a) - U(0) ' ( j 

We further notice that U(h) and U(a) are related to T x and T 2 for the qua- 
sistatic limit J — > 0. In fact, using Eq.(71), we get the relation 

Ujoj ~ U(0) U(h) ~ U(a) = , m 

From Eq.(85) and Eq.(86), we obtain the efficiency for the quasistatic limit 

n = (87) 

This is equal to the Carnot efficiency Recently, Sakaguchi has derived the 
Carnot efficiency by analyzing a model with stochastic boundary conditions [13]. 
His model turns out to be equivalent to that the limit (h — a) — > is taken in 
our model. 



Lastly, we address a few comments about above two examples. In examplel the 
revision term intrinsically contributes to the averaged heat, while the revision 
term has no contribution in example2. This difference is caused by whether 
or not periodic boundary conditions are assumed. As a matter of fact, we can 
generally say that, in one dimensional system, once a steady state is formed, 
the revision term does not affect to the averaged heat under periodic boundary 
conditions. Since the steady condition leads to constant flux by the equation 
of continuity, the averaged value of heat is given by 

< Q H ji /l>' (88) 

This is certainly equal to zero when periodic boundary conditions are as- 
sumed. However, this matter is limited to one dimensional systems. When the 
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dimension is higher than one, 

divJ = 



(89) 



is derived from the equation of continuity. In this case, the circulation u is 
defined by 

u) = rot J, (90) 

which may play an important role in energetics [14]. We expect that the revi- 
sion term contributes to the heat coupled with the circulation. It is a future 
problem to study energetics in high-dimensional non-uniform temperature sys- 
tems. 



6 Summary and Conclusions 

In this paper, we have proposed an energetic interpretation about a stochastic 
process described by a Langevin equation with non-uniform temperature. We 
analyzed two examples as applications of the energetic interpretation. As the 
first example, the Thomson coefficient for non-interacting particles was proved 
to be 1/2. The second example concerns the Buttiker's model. We calculated 
the work efficiency and showed that the Carnot efficiency is realized under the 
quasistatic condition. 
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